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Abstract 

Hydrophobic interactions provide driving forces for protein folding, mem- 
brane formation, and oil-water separation. Motivated by information theory, 
the poorly understood nonpolar solute interactions in water are investigated. 
A simple heuristic model of hydrophobic effects in terms of density fluctu- 
ations is developed. This model accounts quantitatively for the central hy- 
drophobic phenomena of cavity formation and association of inert gas solutes; 
it therefore clarifies the underlying physics of hydrophobic effects and permits 
important applications to conformational equilibria of nonpolar solutes and 
hydrophobic residues in biopolymers. 



^present address: Department of Chemical Engineering, University of Delaware, Newark, DE 
19716, USA 

■^Department of Pharmaceutical Chemistry, UCSF, San Francisco, CA 94143; and NASA - Ames 
Research Center, MS-239-4, Moffett Field, CA 94035-1000, USA 

*To whom correspondence should be addressed. LA-UR-95-3486 



Hummer, Garde, Garcia, Pohorille & Pratt: Hydrophobic interactions 



Page 2 



Protein folding and formation of lipid membranes and micelles are driven by hydrophobic 
interactions. Understanding the nature of these interactions of nonpolar solutes in water and 
quantitatively describing them has been the topic of numerous experimental and theoretical 
studies (|T|). Despite these efforts, a general, quantitative, molecular scale theory is not 
available. Motivated by information theory, we study molecular scale hydrophobic effects 
by exploring the thermodynamics of hydration and association of model hard core solutes 
(0) where otherwise the only reliable calculational tools are elaborate computer simulations. 
As we will show, a simple theory based upon local density fluctuations proves sufficient to 
describe qualitatively and quantitatively the central hydrophobic phenomena. The theory 
can be applied easily to solutes of arbitrary shape and size, and permits study of structural 
equilibria of nonpolar molecules and residues. 

Among the previous theories of hydrophobic hydration, the scaled particle model 
initiated an important line of development. It utilizes asymptotic results in the small and 
large solute limits where direct information on solvent correlations is not required. In the 
intermediate region and for phenomena such as solvent separated hydrophobic interactions, 
(|,^ molecular correlations are important. Stillinger (0) and others (H) partly filled these 
gaps, but much remains missing in our knowledge of hydrophobic effects. Our understanding 
of hydrophobic interactions is particularly meager. The available solubility or "hydrophobic- 
ity" models have not encompassed the molecular structural features that give rise to stable 
free energy minima with water separating hydrophobic groups — the so-called solvent sepa- 
rated hydrophobic interactions. Thus, the hydrophobicity models have not been justified by 
reproducing these interactions and to this extent must be regarded as of unproven validity 
for conformational analysis of biomolecular structure in solution. 

Chemical potentials of nonpolar solutes in water are the quantities of ultimate interest 
in the theory of hydrophobic hydration. We shall mainly be concerned with the excess 
chemical potential Z^^^ of a hard core solute in water. Statistical mechanics relates /S^^^ to 
the probability po of finding an empty volume, or cavity, of given size and shape in water at 
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equilibrium, Z^^^ = — /c^Tlnpo- However, po becomes exceedingly small for larger cavities, 
and direct calculations of po (e-g-; by computer simulations using test particle insertion (j^)) 
become impractical. Our goal here is to provide a model to access this region of cavity 
sizes based upon accessible, even if indirect, information. The model we seek should not 
only extend the range of sizes but also apply to non-spherical cavity shapes. Complete 
information on the binomial moments of the occupancy of a volume identical to the cavity 
of interest would permit a construction of po based upon the well known expression (||,||,^ 
that exactly relates po to the fluctuations in the number n of solvent molecules in the cavity 
region: 

— £m)'(^)- w 

This formula is useful in guiding the analysis of theories (|^ but becomes impractical for 
larger cavities because of its sensitivity to missing information. In those instances, higher 
binomial moments cannot be neglected and that information is unavailable. We require an 
approach that is less sensitive to missing information while exploiting information readily 
available for cavities of all sizes and shapes. These are the low order moments of the number 
n of solvent molecules in the cavity volume Al/, 

(n) = pAr, (2) 
(n(n-l)) =p^ I dv I dv'g{\v-v'\), (3) 

where g{r) is the water oxygen pair correlation function and p is the water density. Infor- 
mation theory (|T0|) provides a paradigm for exploiting such quantities. We consider the set 
of probabilities p„ of finding exactly n solvent molecules in the cavity volume, Y^^=oPn = 1- 
Observation of no solvent molecules in the cavity region is then just one of the elemen- 
tary events, and po is just one of the desired probabilities. We formulate models of the 
distribution p„ that satisfy known moments, Eqs. ^ and ^. Minimizing an information mea- 
sure permits the inference of probabilities of maximum likelihood in a sampling experiment 
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that satisfies the moment information ([TD|). However, the accurately parabohc character of 
lnp„ for cavity volumes of small molecular size, observed in computer simulations of wa- 
ter (0), suggests an even simpler model: We adopt the form pn = exp(Ao + Xin + A2^^) 
with "Lagrange multipliers" Ao,Ai,A2 to be determined by the three moment conditions 
J2'?^=oPn = I5 J2'i^=o^Pn = {n), and J2'?^=o^'^Pn = (^^)- From this we extract the po that 
provides the desired thermodynamic result. 

Figure |l] shows the chemical potential for a spherical cavity in water calculated from 
the model as well as from direct cavity statistics. We find excellent agreement in the range 
considered. This simple theory based entirely on the pair correlation function of water 
is capable of reproducing the thermodynamics of cavity formation in the region that is 
accessible to direct computer simulations (|^,p!^-[l^. 

So far, we considered the hydrophobic hydration of a single cavity. We now proceed 
to study hydrophobic interactions. The free energy associated with bringing together two 
hydrophobic solutes corresponds to a distance dependent potential of mean force (PMF). The 
developed methodology provides us with the means of calculating this PMF. Given the first 
and second moment of the particle number distribution, an approximate chemical potential 
Z^^^(r) can be determined for forming a cavity made of two spheres with given radius R and 
distance r using the model. The PMF is then defined as W{r) = /S^^^{r) — limg^oo ^'^^{s). 

In the following, we study the association of two cavities of methane size in water. For the 
radius of the spheres, we have chosen a value of R = 0.33 nm. This radius corresponds to the 
smallest distance where methane- water pair correlations reach 1.0 in commonly employed 
models ([T5|). The cavity PMF is shown in Fig. |^. As a reference, the cavity potential 
produced by the molecular dynamics simulation of Smith and Haymet (|T6D is included. We 
obtained this by subtracting the bare potential from the PMF of methane association. Again, 
we find qualitative and quantitative agreement between our simple model and elaborate 
computer simulations. Our cavity PMF shows a strongly favored region with overlapping 
cavities, separated by a distinct barrier from a solvent separated, stable minimum at about 
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0.7 nm distance, in agreement with the computer simulations for methane association. In 
addition, we observe a very shallow third minimum at 1.1 nm distance. 

As a last example, we study the torsional equilibrium of n-butane. Figure ^ shows the 
cavity PMF as a function of the torsional angle 0, which is compared to explicit computer 
simulation results of Beglov and Roux (|l3)- In complete agreement with the computer 
simulations, we find that the more compact cis (0 = 0) and gauche structures (0 = vr/S) 
are favored over the extended trans conformation (0 = vr) by about 1.8 kJ mol~^ and 
0.7 kJ mol~\ respectively. 

When the cavities of interest are much larger than the size of the solvent molecules, 
other physical considerations intrude, as is well known in the scaled particle developments 
(0)- We expect this two moment model to become unsufficient as a sole description in those 
cases. Nevertheless, a broad range of problems of biological interest, such as interactions 
of ligands with binding sites, effects of point mutations on protein-solvent interactions, and 
conformational equilibria of side chains, is within the range of applicability of this model. 

The present theory permits a sound, quantitative understanding of hydrophobic hydra- 
tion and association processes. This view is based upon (i) the simplicity of the present 
theory involving only particle number fluctuations, (ii) the small amount of input informa- 
tion required [g{r)\ that is experimentally accessible, (iii) the computational ease of eval- 
uating the theory permitting application to realistically large and complex shaped solutes, 
(iv) the unified and consistent treatment of hydrophobic hydration and interaction phenom- 
ena, (v) the clear connection of information theory to basic statistical thermodynamics, and 
(vi) the qualitative and quantitative accuracy of the theoretical predictions in a range that 
is otherwise inaccessible to direct theoretical treatment. 
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FIGURES 

FIG. 1. Hydrophobic hydration: Comparison of the chemical potential calculated from simula- 
tion (symbols) and the fluctuation model (line) for spherical cavities with radius i?, the distance 
of closest approach of a water oxygen to the solute. The fluctuation integral Eq. ^ required by the 
information theory was reduced to one-dimensional integration (|l8|). The simulation results were 
gathered from test-particle insertion, where 8000 configurations (separated by 50 passes each) of a 
Monte Carlo simulation of SPC water (^) were used. The same simulation was used to calculate 
the water-oxygen pair correlation g{r). 

FIG. 2. Hydrophobic interaction: PMF of cavity association, r is the distance of two spherical 
cavities with radius R = 0.33 nm. The result of the fluctuation model (solid line) is compared 
with the cavity PMF of Smith and Haymet (dotted line; from Fig. 4 of Ref. ([l^). The inserts 
illustrate the size of the excluded volume. At the solvent-separated PMF minimum (0.7 nm), a 
water molecule barely fits between the two hard-sphere solutes. 

FIG. 3. Torsional PMF of n-butane. Butane was modeled as four spheres with radius 
R = 0.33 nm, bond length 0.153 nm and tetrahedral bond angle. The result of the fluctua- 
tion model (solid line) is compared with the cavity PMF of Beglov and Roux (dotted line; from 
Fig. 8 of Ref. ([TtI)). The cis, gauche, and trans rotational states are shown as inserts. 
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